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AT SUPERSONIC SPEEDS 
(AN EXTENSION) 

By John C. Evvard 


SUMMARY 

A derivation is presented of a point-source method based on 
the linearized theory for obtaining the pressure coefficient on 
thin wings at supersonic speeds. The method is applied to calcu- 
late the lift distribution of a thin flat -plate wing having a 
straight leading edge and an arbitrarily curved wing tip whose 
elements are swept behind the Mach angle. A qualitative basis for 
obtaining solutions that satisfy the Kutta-Joukowski condition is 
included . 

The analysis is continued to formulate the velocity potential 
for regions influenced by so-called subsonic trailing edges (that 
is, edges where the component of the free-stream velocity normal 
to the trailing edge is subsonio) . The derivations include the 
solution that satisfies the Kutta-Joukoweki condition or any of 
the Infinite number of transition solutions involving a disconti- 
nuity in the cross velocity (sidewash) in the wake of a subsonic 
trailing edge. As an example, the two perturbation-velocity 
components in the plane of a trapezoidal wing and the upwash over 
the wing tip are evaluated. By means of series expansions, the 
mathematical nature of the lift distribution in regions influenced 
by arbitrary subsonic trailing edges is indicated when the Kutta- 
Joukcwski or some other condition is applied to make the solution 
unique . 


- INTRODUCTION 

The problem of determining the lift distribution of arbitrary 
thin wings at supersonic speeds has been partly solved by use of 
the linearized theory in references 1 and 2. In reference 1, the 
point-source method of reference 3 is extended to include the 
effective sources generated by the slopes of the streamlines in 
the region between the wing boundaries and the foremost Mach waves. 
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The velocity potential and hence the aerodynamic coefficients in the 
vicinity of wing tips at supersonic speeds can then be calculated. 

An explicit solution for the slopes of the streamlines in the 
region included between the wing boundaries and the foremost Maoh 
waves is presented in reference 2. This solution allows the methods 
of reference 1 to be extended to include the calculation of the 
velocity potential at points on the wing influenced by interacting 
external flow fields. 

In subsonic aerodynamic theory, an infinite number of solutions 
will satisfy the boundary conditions for a given wing. Some 
criterion is then required to choose from the manifold of solutions 
the one that is experienced in practice. A criterion that yields 
theoretical results in approximate agreement with experiment is the 
Kutta-Joukowski condition. This condition requires that the stream- 
lines leave the trailing edge of the airfoil smoothly (that is, that 
the three components of the perturbation velocity are continuous 
across the trailing-edge boundary) . 

When the flow is supersonic, the Kutta-Joukowski condition 
need not be fulfilled at the trailing edge of the airfoil. Local- 
ized compression and expansion waves equalize the pressures from 
the bottom and top surfaces behind the trailing edge. Nevertheless, 
because of the close similarity of supersonic and subsonic flow 
when the component of the free-stream velocity normal to the 
trailing edge is subsonic (the so-called subsonic trailing edge), 
a strong opinion among aerodynami cists is that the Kutta-Joukowski 
condition will also approximately hold at supersonic speeds. 

As indicated in reference 2, the formulations of reference 1 
may be extended to include solutions that satisfy the Kutta- 
Joukowski condition. The extension shows that in supersonic-flow 
theory, as in subsonic-flow theory, an infinity of solutions can 
exist for a given wing boundary. Some additional criterion is 
then required. For subsonic leading edges, the added boundary 
condition applied in reference 1 is that the perturbation-velocity 
components shall be continuous in the region between the wing 
boundary and the foremost Mach waves. The solutions so obtained do 
not satisfy the Kutta-Joukowski condition for regions influenced 
by subsonic trailing edges. If a discontinuity in the oross- 
velocity component (sidewash in the plane of the wing) is per- 
mitted (corresponding to the shedding of a vortex sheet), solu- 
tions may be obtained that satisfy the Kutta-Joukowski condition 
or any transition solution depending upon the assumed strength of 
the discontinuity in the cross velocity. The particular 
criterion to be used can be determined only by experiment. 
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A derivation based on linearized theory is presented in the 
first part of this report to show that the theoretical pressure dis- 
tribution on thin wings at supersonio speeds may be obtained from, a 
line integral and a surface integral. (A line source was applied in 
reference 4 to obtain the pressure distribution at zero lift.) For 
thin flat-plate wings or for the lift distribution on finite- 
thickness wings, only line integrals are required. As an example, 
the lift distribution is evaluated for a wing having a straight 
leading edge intersected by an arbitrarily shaped curve whose 
boundary-line elements are swept behind the Mach line from the 
wing-tip and leading-edge intersection. The second part of the 
report formulates the expressions for the velocity potential in 
regions influenced by subsonic trailing edges. 

Examples of these methods are presented. The solution satis- 
fying the Kutta-Joukowski condition is derived for a swept trap- 
ezoidal wing. Sidewash and upwash velocities in the vicinity of 
the wing tip are also calculated. As a further application, by 
* means of series expansions, the mathematical form of the lift dis- 

tribution associated with the wake behind a general subsonic 
trailing edge is indicated. When this result is combined with the 
line integrals in the first part of the report, expressions are 
obtained for the lift distribution of arbitrary-plan-boundary 
single wing tips at supersonic speeds. The Kutta-Joukowski and 
other conditions are applied to formulate feasible solutions for 
the various wing areas. The research was conducted at the KACA 
Cleveland laboratory during August and September 1947. 


ANALYSIS 


Point-Source Method for Evaluating 
Pressure Coefficient 


A convenient approach toward evaluating the pressure coeffi 
cient on thin finite wings at supersonic speeds is to establish 
the velocity potential. A point-source integral expression for 
the velocity potential at any point on thin wings at supersonic 
speeds is given in references 1 and 3 as 


q> 


u 

Jt 


T 


qdfdq 

'V (x-l ) 2 -P 2 (y-Tj ) 2 


( 1 ) 
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where here 

cp perturbation- velocity potential on top wing surface 

U free-stream velocity considered parallel to x or | axis 

C slope of streamlines near x,y plane measured in y = constant 

planes 

|, r) Cartesian variables of integration, x and y directions, 
respectively 

x,y Cartesian coordinates of point where cp is evaluated 
p = ^M 2 -l 

M free-stream Mach number 

(A complete list of symbols used in this report is presented in 
appendix A.) In equation (l), the integration of a is conducted 
in the x,y plane for the disturbed flow field S included in the 
forward Mach cone from the point (x,y). For example, the integra- 
tion for the wing shown in figure 1(a) is over the area S^. In ** 

regions influenced by a wing tip, the disturbed flow area between 
the wing boundary and the foremost Mach line must also be included. 

If the velooity potential Is computed at a point — (x+dx,y), as 
in figure l(b), the area enclosed by the curve efgc' is the same 
as the area abdc and the relation of the distances in the denomi- 
nator of the integrand Is the same as was considered at point (x,y) . 

The slopes, however, will have the value a + -g-g dx. In addition, 

the area a'bd'gfe will contribute to the velocity potential. The 
value of cp then becomes 



where f ]_ is the value of | along the leading edge. Integration 
yields 
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where abd represents the path of integration. Subtraction of 
equation (l) from equation (2a) gives 



U 

it 



J(x-|) 2 -3 2 (y-r 1 ) 2 


_ 2 

it 

■J abd 

where £ and CT are evaluated in the line integral as functions 
of q along the line abd of' figure l(a). 

The x component of the perturbation velocity may then be 
computed as the surface integral of the rate of change of the stream- 
line slopes with respect to £ plus a line integral along the 
foremost boundary of the disturbed flow field included in the 
forward Mach cone from the point (x,y). The surface integral, of 
course, vanishes in regions where O is constant. (Evaluation of 

the surface integral along lines of discontinuous ^ results in 
additional line integrals . ) 

For some problems, the integrations of equation (3) may be 
more easily evaluated in a set of oblique coordinates (u,v) whose 
axes lie parallel to the Mach lines. The transformation equations 
are 
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u - 

V 

-§«+ n) 

i = |(v+u) 


- 5 (t - u) 

u v * 2p-(3c - Py ) 

v w 

- §(x + py) 

X - f(v w -Hlw) 

J 

3 ^ V v-«w) 


( 4 ) 


Inasmuch, as the elemental area in the (u,v) coordinate system 
is -^Idudv, equation (3) may be written as 


M‘ 


U 

2p« 


U 

2 Pit 


[f fe - 

Js ^(uv-u) (v w -v) 

g (dv-du) 
abd ^(uv-u)(v w -y) 


(3a) 


where (u w , v w ) corresponds to the point (x,y) and the evaluation 
of the line integral' Is again conducted along the foremost line of 
disturbances included in the forward Mach cone. 

When only one wing tip is included in the forward Mach cone 
from the point, the surface integral over the area Sp (fig. 2(a)) 
generated by the upwash over the wing tip may be replaced by an 
equivalent integration over a portion of the wing. This simplifi- 
cation, which follows from the derivations of reference 1, 
eliminates the necessity of evaluating the slopes of the streamlines 
in the external flow field Sq. The velocity potential at point 
(uw,v w ) cm the wing of figure 2(a) Is obtained from reference 1 as 

g rpdudv 
tyOv-uKvv-v) 

(gj»-tgm ) dudv 

< 5 > 

2 ^(iv-u) (vw-v) 
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where <jb and a ij are the bottom- and top-surface wing slopes, 
respectively. 


Because evaluation of cp from equation (5) requires the applica- 
tion of separate sets of wing slopes in the respective regions 
and S V) 2) the calculation of rg will result in a line integral 

along the dividing boundary. The derivation presented in appendix B 
then yields equation (6) . 



U 

2 Pit 


®w,l 


'Sc T dcr T \ 
vST + w dudir 

^(u^-u) (v w -v) 





S(Cfg+Gji) 

Su + Sv j, j, 
— dudv 

(l^-tt) (v w -v) 


_y_ 

2 pit 


r- 

^ab 


qiji(dv-du) 
/\/(uw-u) (v w -v) 


r- 

II 

2 Pit 

w»bod 


(Ojj-KJrji) (dv-du) 

2 hj (u w -u) (v w -v) 


IT 

(VV 

• 

du 2 ( V 
L " dv w J 

dv 

2 Pit 

bd » 2 V 0v u > ( v w-v) 



( 6 ) 


du 2 (T w ) 

where — 

dv w 

at v = v w . 


\ 

is the derivative of the wing-boundary curve u 


U2(v) 


In order to obtain the lift distribution as a function of angle 
of attack a, only the solution for a thin flat plate is required. 
Because ag = -a^i = a, only the two line integrals along the curves 
ab and bd remain. The portion of the equation (6) that yields the 
lift distribution is then 
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r 


U 

2P* 


Jab 

n- 


u 

2 Pit 


bd 


-a (dv-du) 
tfiuv-u) (v w -v) 

- a ( 1 - 5 ^V T 

/ V( u v" u ) (▼*-▼) 


(7) 


If the equations of the wing boundaries (fig. 2(b)) on the two 
branches separated by the origin are denoted by 


= vi(u) or u = u^(v) 


and 


v = v 2 (u) or u = u 2 (v) 
equation (7) may be partly evaluated to yield 


s?- 


r 


TJa. 
2 Pit 


(dv-du) 




r 


ab 


^(Uy-u) (V W -V) 


Ua 

Pit 


Cl- 


du 2 (v w ) 


dv. 


w 


Vy-vi(u 2 (v v )) 

Uw“W2(v w ) 


(7a) 


where v^ Is a function of u 2 (v w ) . The expression for ■^ > 

represented by equation (7a) Is applied in the second part of this 
report to derive solutions for the lift distribution that satisfy 
the Kutta-Joukowski condition for general-plan-form single-wing- 
tip boundaries. 

As an illustration of equation (7a ) , the pressure coefficient- 
on the top wing surface as given by the relation 


C P ' 


( 8 ) 


will be calculated for the wing shown in figure 3. 
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Inasmuch as the leading edge and the wing tip satisfy the equa- 
tions v = - k]_u and u = u 2 (v) , the integral of equation (7a) 
becomes 


Ua 
2 Pit 


. u 2(v w ) 


(l+k]_)du 


u w 


V(uv-u) (v w +k!u) 


Ua 


(* 1 + 1 ) 


tan' 


-1 


ki[u v -u 2 (7v)] 

Vv+kiuaCr-w) 


and 


Ua 

P« 


1- 


du 2 (v w ) 


dv. 


w 


T w -vi(u 2 (v w )) ua 


\ «w-^2( v w) 

The value of Cp is then given as 


“ pn * 1 “ 


d u-2( v w) 


dv. 


v j 


Vy+kj^CVy) 

uv-u 2 (v w ) 


( 9 ) 


( 10 ) 


2a 

Pn 


r ^(^w)] 

Ivv+kiu2(v w ) ki+1 

,|ki[uw-u 2 (v w )] 

[L dT » J 

^ iV-u 2 (v v ) 

^ v w +kiu 2 (v w ) 


( 11 ) 

For comparison, equation (11) is also derived in appendix C as equa- 
tion (C4) by the method of reference 1. 


Equation (ll) shows that the pressure coefficient assumes the 
Ackeret value for a swept wing along the innermost Mach line from 
the origin (v w = 0) and approaches infinity along the subsonic leading 


edge of the wing tip. The portion of the solution that gives the 
infinity (that is, equation (10) ) is generated by the line integral 
along the Mach line bd in equation (7). This quantity is zero when 

the sweep on the wing tip is -r or — = 1. A change in sign 

du 2 . 

occurs when vr - > 1 (regions influenced by subsonic trailing edges) 


dv. 


w 


in equation (10), which leads to a reversal of sign in the pressure 
coefficient over a portion of the top wing surface. This pressure 
reversal leads to lost lift and, because of trailing-edge suction, to 
decreased lift-drag ratio as compared with the solution that satisfies 
the Kutta-Joukowski condition. 


The pressure coefficient that satisfies the Kutta-Joukowski con- 
dition along the subsonic trailing edge must match the solution given 
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in equation (ll) along the oonatant-v w line evaluated by the 
dU2 

relation — = 1. The pressure coefficient must furthermore he 

zero along the subsonic trailing edge. A function that meets both 
of these requirements is 


2a ( k l +1 ) ^ n Wki[u w -U2(vv)] 
fjkl ** ( v w +kiu 2 (v w ) 


( 12 ) 


The wake from the subsonic trailing edge alters the lift distribu- 
tion over a portion of the wing in such a manner that the discon- 
tinuous tern in equation (ll) is canceled. In the second part of 
this report, equation (12) is shown to be the correct solution for 
the pressure coefficient in regions influenced by subsonic trailing 
edges of a wing whose leading edge is straight (fig. 3) when the 
Kutta-Joukowski condition is to apply. 


Formulation of Velocity Potential Including Effects 
Associated with Wake behind 
Subsonio Trailing Edge 

In order to clarify the present discussion, a brief review of 
the concepts presented in reference 1 is included. The calculation 
of the velocity potential at some point on the surfaoe of thin wings 
at supersonic speeds requires a knowledge of the streamline slopes 
( a of equation (l) ) included in the forward Mach cone from the 
point. The slopes of the streamlines over the area S^i+g) 

(fig. 2(a)) are Just the wing-section slopes. The slopes of the 
streamlines in the region Sj) must then be calculated, either 
implicitly or explicitly. 

The slopes of the streamlines in the region Sp may be evaluated 
by assuming that a thin impermeable diaphragm coincides with a stream 
sheet in the plane of the wing. This assumption does not alter the 
flow over the wing. The diaphragm may then be regarded as an exten- 
sion of the wing to isolate the top and bottom surfaces. 

The coincidence of the diaphragm with a stream sheet requires 
that no pressure difference exist between its top and bottom 
surfaces. This condition may be satisfied by requiring either that 
the surface velocity potential be continuous in passing from the top 
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to the bottom surface of the diaphragm, as in reference 1, or that 
the component of the perturbation velocity along the free-stream. 
direction (the z component) be continuous. 

If the x component of the perturbation velocity (obtained by 
differentiating the velocity potential) is equated on the top and 
bottom surfaces of the diaphragm, the defining equation for the 
slope X of the stream sheet S B (fig. 4(a)) is obtained as 


5 

5E 


rr 


M£,n)dSdT| 


%> 


ty(x-t) 2 _( 32 (y_ri )2 


a 

= SE 


(c B -a T )d|dTj 




2 


(13) 


Integration of equation (13) yields 




X(|,tl)d|dTi 


rr 


S D 


V(*-£) 2 -0 2 (y-n) 2 


— 2 — 


d|dT) 




£) 2 -3 2 (y-T}) 2 


|f(y) 


(14) 


If the function f(y) that results from the partial integration is 
zero, the equation derived in reference 1 for defining X is 
obtained . 


A calculation of the velocity potential on the top surface of 
the diaphragm for any value of the function f(y) may be made. 
Again with reference to figure 4(a), the velocity potential becomes 


f 

u 


a T d|dTi 


rt 

'■A. 


|) 2 -P 2 (y-n ) 2 


r 




u 

■ 

Xd|dTj 

(15) 

Jt 

J. 

Jsj) 

^(x-D 2 -p 2 (y-n) 2 



Substitution of equation (14) into equation (15) yields 



(a B -K7 T )d|dTi 

2 ^(x-S) 2 -P 2 (y-Ti) 2 


+ f(y) 


(16) 
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A similar calculation of the potential on the bottom surface of the 
diaphragm changes only the sign affixed to f (y) . The function f(y) 
thus represents one-half the difference in potential between the top 
and bottom surfaces of the diaphragm. In the case of thin plate 
wings, = 0 and f(y) is the potential on the top surface of 

the diaphragm . 

Despite a discontinuity in the velocity potential across the 
diaphragm, the x and z components of the perturbation velocity 
are continuous. This apparently contradictory state of affairs may 
be rationalized by the. presence of a vortex sheet in the plane of 
the wing. The strength of the vortex sheet may be chosen (by 
controlling f(y)) to yield solutions for the velooity potential 
that satisfy the Kutta-Joukovski condition for subsonic trailing 
edges . 


The function f(y) is related to the strength of the vortex 
sheet. If no wing nor material objeot is intercepted (for the 
values of y under consideration) to generate a wake, the strength 
of f (y) remains unaltered for all values of x. 

In the equations that follow, two simplifications are introduced: 
(l) all the derivations are conducted in the oblique coordinate 
system, and (2) the derivations consider only thin flat-plate wings.- 
Neither simplification decreases the generality of the solution. 

The results for the actual wing at angle of attack may be evaluated 
as the sum of the results at zero angle of attack and the solution 
at angle of attack for the thin flat" plate of the same plan form. 

The external flow field may be divided into two parts (as in 
fig. 4(b)), separated by the value of y that- denotes the Junction 
between the leading and trailing edge. The foremost Mach line 
originating on the leading edge generally represents a line of - 
infinitesimal disturbance along which f (y) may be set equal to 
zero. The function f(y) will remain zero for values of x that 
do not intercept the wing. The velocity potential (in the 
z = 0 plane) is thus zero in the field Sp^g* The field Sp^p, on 
the other hand, lies in the wake of the wing. In this region, f (y) 
will generally be other than zero (corresponding to the generation 
of a vortex sheet) and may be evaluated along the wing trailing 
edge. 


If the region Sp p temporarily considered as a portion of 
the wing, the method o£ reference 1 may be applied to evaluate the 
potential on the top surface of the diaphragm Sp ]_ or on the wing 
surface. The virtual wing tip is then the y = constant line 
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denoting the junction between the leading and trailing edges. Inas- 
much as f (y) is zero in the region Sj^ 2 , the effects of this 
disturbed flow field may be evaluated by the methods of reference 1$ 
equation (5) applied to the virtual wing boundary will result. When 
the flow fields are subdivided and renumbered, as on figure 4(c), 
the contribution to the velocity potential of the disturbed 
fields Sd( 2-*4+3) and Sw,2 ls given by equation (5) as 


Mrt 




Sw, 2 


(a B +CT T )dudv tj 

2 /y (u^-u) (v w -v) m 




®D,3 


^3-^g)dudv 
2 /\J (Uy.-u) (v w -v) 


= 0 


because a B = -a T . Only the slopes of the streamlines in the 
regions S^i and Sp need then be considered. The velocity 
potential thus becomes 


rr 


f(y) = - 


j u_ 

Mit 


-ctdudv 


U_ 

Mn 


rr- 

s- 


(u D -u) (v D -v) 


Xdudv 


>D,1 


l\j( u D -u) (v D -v) 


(17) 


The origin of the coordinates for the wing of figure 4(d) is 
placed at the junction of the subsonic leading and trailing edges. 
The wing boundaries are defined on the two sides of the origin by 
the two sets of equations: 

v = V]_(u) or u = U]_(v) 

v = v 2 (u) or u = u 2 (v) 

du 2 

(It will also be implicitly assumed that > 1.) Insertion of 

the limits of integration into equation (17) yields 


f(y) 





— V 2 (u) 

f'Vjj 

fVD-UpA U 

du 


-adv 

+ 

Xdv 

\ M / Mrt 







A) u D -u 
v D » 

•o 

V 1 (u) ^ J 

T*(U) 


( 18 ) 
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w- 


By use of Abel's solution (reference 2, 5, or 6), equation (18) may 
be inverted to yield 



The velocity potential at points on the wing influenced by the 
subsonic trailing edge (fig. 4(e)) is 



The second member of equation (19) may replace the first member 
along lines of constant vp (or v w ) that extend across the wing. 
Equation (20) then becomes 



Equation (2l) gives the velocity potential on the wing in regions 
influenced by subsonic- trailing edges . Inasmuch as 

f(y) = f may he arbitrarily chosen, it is apparent that an 

infinite n umb er of solutions can satisfy the boundary conditions. 
If f = 0, equation (21) reduces to the result that would be 
obtained by the methods of reference 1. On the other hand, if the 
Kutta-Joukcweki condition is to apply, the function f may be so 
chosen that - the flow leaves the subsonic trailing edge smoothly. 


After the function f has been chosen, the second member of 
equation (19) may be evaluated. The slopes of the streamlines A may 
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then, he calculated by use of Abel's equation. These slopes in 
general consist of two parts, X = X a +A b , which satisfy the two 
equations 


r v 2 (ti) 


-adv 


rv- 


D X-^dv 


^v-l(u) 


Jv 2 ( u) ^ 


= 0 


( 22 ) 


and 


r VD \dv 


v 2< 


u)f^ 


M d 
U 3u 


^u 


f (y)du D 


v D f ^ 


(23) 


Equation (22) has been evaluated in reference 2 to give 


2a. 

A 

v 2 (u)-v 1 (u) n 

v 2 (u)-vi(u) 

n 

n 

v-v 2 (u) ^ I 

v-v 2 (u) 


(22a) 


Equation (23) may be solved by Abel's equation (reference 2, 5, or 6) 
to give 





f~U /_ 

M d j 

r . 

s 

OUL^ 

*1 

IS—1 


^Vp(u) 


JEEI 


du D 


dv T 


(23a) 


In this maimer, the slopes of the streamlines in the region Sj, ^ 
of figure 4(b) may be evaluated. If the region Sp ± is treated 
as a portion of the wing, the slopes of the streamlines in the 
region Sj ^ 2 te evaluated from the equations of reference 2. 


Examples of Solutions That Satisfy 
Kutta-Joukowski Condition 

Swept trapezoidal wing, - As a simple example, the theory of 
the trapezoidal wing shown in figure 5 will be outlined. The leading 
and trailing edges of the tip are defined by the equations 
v = V]_(u) = - k]_u and v = v 2 u = k 2 u, where k^> 0 and 0 < k 2 < 1. 
Because the wing boundaries are straight lines, conical flow will be 
assumed . 
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Because conical flow has been assumed in the region influenced 
by the wing tip, the cross velocity ^ must be either consent or 
a function of y/x. Because f is a function only of y, ^ ^ 

is constant in the region Sd # i. The function f(y) is then of the 
form aQ + a^y where aQ and a^_ are constants. Inasmuch as f(y) 
must be zero along the line y = 0, ao is zero. Equation (19) then 
becomes 


^ 2 u 




-adv 




^ dv 


M 5 

= " U <5u 


a 1 (v D -u D )du D 2a x 

— : — — 'p'S 

iv D M /yu-uu 


Lk lU ^ Jk 2 u ^ 

Equation (21) then yields the velocity potential as 


(24) 


cp = - 



r v w/ k 2 

r u w 

r v w 

2a^ 

- U 

du 

-adv 

1ST 

r w K- a “J 

v„/lc 2 

-*!» / il T ’'- T 


(25) 


The two integrations indicated in equation (25) may be evaluated by 
formulas 111 and 113 of reference 7 to give 

C Y v/ k 2 .1 

^«- v w , 2a l 

■ . ■ nia = . . 


2a l 

TtT 


V v 


+ (v w -u w ) tan' 


Mrt 


il-k 2 

kc 


^’ v w( ]c 2 u w“ T w) 


^ A N (1 - k2) 

\ k 2 u w -v w 


(26) 


and 


CPb = - 


U_ 
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'U, 


l w 


du 


v w /* 2 r^. 
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-k^u 


-ad’v _ 2Ua 
— = Mit 


V 


'v 


* ^v w ( k^-Vv) 


+ tan -l 




ki(k 2 u v -v v ) 


v w (ki+k 2 ) 


(27) 
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The potential <p = q? a + must assume the value aqy along the 
curve v v = k2 u w Equation (25) (from equations (26) and (27)) 


meets this requirement. Moreover, ^ must he equal to aq along 
the curve v w = kgu-jj- if the Kutta-Joukowski condition is to he 
satisfied. Differentiation of equations (26) and (27) yields 




The cross velocity on the wing is then 



(31) 


In a similar manner, 

(32) 


which is in agreement with equation (12) and with the results of 
reference 8. ¥. D. Hayes of North American Aviation, Ino., has 

obtained equation (32) from conical-flow relations. 

Both the x and y components of the perturbation velocity 
are continuous (equations (31) and (32)) across the subsonic trailing 
edge, as required by the Kutta-Joukowski condition. The upwash over 
the wing tip is also continuous. From equation (22a), 


^ _ 0 £ 
dx "5x 
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, 2a 


1 


u D (k 1 +k 2 ) 


vp-k 2 u D 

From equations (24) and (23a ) , 

X - d 
Aa T5T 57 


- tan' 




u D (k 1 +k 2 ) 


VD-k 2 uj) 


(33) 


rv 


Jk 2 u 


u-vp 

^-v D 


dVD 


Evaluation of X a and replacement of u and v by u D and vp, 
respectively, gives 


% , , 2a i k l +k 2 i 

X a ( u I»TD) = " ~ IfT^ ] 


1 . up(l-k 2 ) 


^v D - 


k 2 uc 


+ log 




D 


^v D -k 2 UD + ^up(l-k 2 ) 


+ 2 " 


V D _U D 


2 *|vp-k 2 up ^A|v D -k 2 u D + /^up(l-k 2 ) 
Addition of X a and X b yields 


(34) 


. 2a 

A — - ' 

it 


tan 


■iij 


(k^+k 2 )up ^|k 1 +k 2 

v D -k 2 u D 


A K l +K 2 

\ l“ k 2 


log 




^vp-k 2 u D + up(l-k 2 ) 


(35^ 


Equation (35) is valid only in the region Sp ; i of figure 5. As 
v approaches k 2 u, X from equation (35) approaches -a, as 
required by the Kutta-Joukovski condition. The solution then satis- 
fies the Kutta-Joukovski condition in all the perturbation-velocity 
components even though only one of the components was applied to 
evaluate the constant a]_. 


The x an d y components of the perturbation velocity and the 
slopes of the streamlines in the region Sp^i are plotted as a 
function of Py/x in figures 6 to 8 for a ving on vhich k x = 2 
and k 2 = 0.5. The x and y components (figs. 6 and 7) of the 

“ /x is changed from -1 (inner- 
Both the x and y com- 

_ i< ilz o. 

3 x 


perturbation velocity decrease as Py/a 
most Mach line) to -l/3 (the ving tip). 

ponents remain constant in the region 


The streamline 
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slopes (fig. 8) are equal to -a on the wing. The slopes change 
from -a to +» in the region - ^ ^ ^ ^ 0. Similar values were 
found by Lager strom in an investigation for the Douglas Aircraft Co. 


Generalized wing tip. - For this example, a function f(y) that 
can he expanded in a power series may he chosen for the region Sp^i 
of figure 4(h). The function f(y) may then he expressed as 


f(y) 




= a 0 + 


a l( T p“ u p) a 2^ v D~ u D^‘ 


M 


M 


(36) 


Also, from formula 750 of reference 7, 

1 1 _ H, 1 /TD~ u D \ 1X5 fyD-upV 

2 V u-Vp / 2X4 ^ u-vp/ 


u-up ^Ju-vp 


1X5X5 

2X4X6 


pCSf + . . . 

\ u-v Dy / J 


(37) 


The In tegrand in the second memher of equation (19) then "becomes 


fM L_ 

\ U-Up 


a 0 + 


a l _ a Q 

M. 2(u-v D ) 


J (vp-up) 


a 2 a l 1X3 I /_ ..^2 

aj_ 1X3 &q 


a 3 


a 2 


M 3 2M 2 (u-v d ) M 2x4(u-vp)' 


ap 1x3x5 
2x4x6 (u-Vp) 3 J 


(vp-up ) 3 + . . . L 


(38) 
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Integration of equation (38) with respect to uj> yields 

1 

&o(u-v D ) 2 ■ 


ru 1 

f (y)du D 


V D 


fju-u 


■D 


f2 a l 1x5 a 0 

m 2 " 2M(u-v n ) + 2X4 


a 0 

3 

(u-vd) 2 

2 (u-V£r) 

2 


(u-v D ) 


[*3 

V " 


a 2 


1X3 &! 


(u--qj) 2 

3 

7 

1X3X5 

a 0 

(u-v D ) 2 

2X4X6 (u-Td)^J 

4 


(39) 


Equation (39) may "be simplified by collecting the coefficients of 
the constants ag, a^, ag, etc.: 


ru 


V D 


f (y)du D h ( 1 

®o (u-v D ) 2 ( 1 + 

^U-V D \ 

(u-t D ) 5 (| ♦ 


. 1X5 1X5X5 

2X2 2X4X3 2x4x6x4 * * 


fa 

M 


1 1X5 1X5X5 

2X3 2X4X4 2X4X6X5 * 


+ (u-v D ) 2 


[3 2x4 


1X5 1X5X5 

+ 2X4X5 2X4X6X6 + * 


•) 


f3 , 1 , 1X5 , 1X5X5 . \ 

“ m 3 ^ ^ 4 + 2X5 2X4X6 2X4X6X7 * * * ) 


+ • • • 

Now by formulas 750 and 482 of reference 7 , 
r 1 - • - 


x^-^dx _ 1 


0 ^ 


1 1X5 . 1X3X5 

+ *-n'\ i 


( 40 ) 


n 2 (n+1) 2x4 (n+2) T 2x4x6 (n+3) 


1 “' • • 


■n / \ n 
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Equation (40) then becomes 

ru 


1 3 

= a 0 (^-v D ) 2 Xy-| (u-v D ) 2 X 3 ^ 


YD 




2 2*2 


+ (u-t d )2 x 

w- 


2X1 a 3 


_ 22. ( u ^2 v 3X2X1 

3 (u vdJ X ? 5 3 i + • • • • 


2 2 2 


M' 


2 * 2 * 2*2 


CO 

I 

m=0 


2m+ l 

2 r( n + |) 


(I) 


(42) 


Partial differentiation of equation (42) with respect to u yields 
ru 


a 

5T 


-Vi 


ffyJduj) a 0 


D 


ff-^D 


1 

,2 


(u-v D ) 

2 


2*2 


a 3 


M" 


5 

( 2 5X2X1 
5 3 1 
2 * 2*2 


- -=• (u-vj)) 1 - r-^rr + • • • 


CO 

■I 


m=0 


c-D” is (u-vd) 2 — 

*” r (» + 5) 


(43) 


Equation (43) shows that either member of equation (19) is a function 
only of the quantity u-vp. 
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The portion of the velocity potential associated with the 
function f(y) is given from equation (21) as 



where v w has replaced vjj inside the "bracket, because only u and 
u D are variables in the manipulation. From formula 750 of refer- 
ence 7, 



Substituting equations (43) and (45) into (44), with v w replacing 
Vj) and U 2 (v w ) shortened to U2, yields 


88fi 
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Collecting the coefficients of the constants sq, a^, a .2 . . . 

gives 



24 


HACA TN No. 1585 



S*ram a comparison of each infinite series of equation (46a) with the 

series expansion for ( l-x^) 2 (formula 754 of reference 7), the 

following equation may he deduced: 
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(where s Is simply an integration variable) or 

r (n+l)r (|)an(v w -u v )n 

r K)^ 


rtcpa 



u 2- v w 

1%-Vv 


n=0 


3 2n 


ds 


(46b) 


(The integrals of equation (46b) may be expressed as incomplete 
beta functions.) 


Along the boundary U2 = the upper limit of the integra- 

tions becomes unity and the values of the integrals are given by 
formula 482 of reference 7 as 


s 2 * 1 ds _ r 
l\l- s 2 2T(n+l) 


0 
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In this spepial case, then, cp a reduces to 


*c? a 



n=0 


M 11 


(v w -Uv) n 


= itf(y) 


as is required along the curve U 2 =u v . 


The portion of the x component of the perturbation velocity 
associated with f (y) may be obtained by partial differentiation 
of equation (46b) to give 



n=0 


(-l)T( n+ l)r(|)^ n -i 

■ 1N-n W (u S -t») 2 


(47) 


The series of equation (47) is a function only of v w . Alterations 
of the function f (y) for a given plan form cannot change the factor 



/\|«w~ u 2 


When an arbitrary strength of vorticity is allowed, the x component 
of the perturbation velocity may be obtained by addition of equa- 
tions (7a) and (47) : 

dv-du 

^(u^-u) (v w -v) 



XJa 


HD 

^V^2 


00 


<|vn("2) - 





n- 


n=0 


'H) 




(u 2 -v w ) 


(48) 


or 


TJa 

26rt 


^ab 


i 1 - 

(uy-u) (v w -v) \ U W - u 2 


dv-du 


(48a) 
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where g(v w ) represents times the traced portion of equation (48). 

The function g( v w ) depends on the shape of the wing-plan boundaries 
(fig. 4(f)) and on the amount of vorticity in the wake of the sub- 
sonic trailing edge. The quantity 
lines of constant v w . 


1- 


d 


is constant along 


If solutions that satisfy the Kutta-Joukcwski condition are 
desired, the constants of equation (48) may be evaluated. Along the 
subsonic trailing edge, the x component of the perturbation velocity 
must be zero. The integral of equation (48) is zero along the boundary 
because the limits of integration then coincide. The Kutta-Joukowski 
condition can then be satisfied only if 


Al T w-TL( u 2( v v)) = 



n=0 


(-i) n r(n+i)r(|) 

r K) 


^ (u 2 -v v ) n 2 (49) 

M n 


Equation (49) allows all the coefficients an to be chosen. Because 
the first member is finite when u 2 = v w , bq = 0. The evaluation of 

of the other coefficients may be simplified by setting p = |Ju 2 -v w . 
Equation (49) then becomes 


2Ua 

M 


^ w -vi(u 2 (v w )) = 



n=l 


(-i) n r(n+i)r(|) 

r H) 


55. p 2n-l 
M n 


(49a) 


Equation (49a) is a power series in p, whose nth coefficient an 
may be evaluated in terms of the (2n-l)th derivative about p = 0. 
Differentiation of the first member of equation (49a) may be 
accomplished by successive applications of the relation 



2 Aju 2 (v v )— v w 
* ~du 2 (v v ) ~ 

dv w 



(50) 


By application of equation (50) and L'Hospital's rule to equa- 
tion (49a), the coefficients ai and a 2 , for example, have been 


evaluated . 


a l 


-Ua 


\ 


1 


dV]_ du 2 

du 2 
dv w 


1 


lim u 2 - 4 v w — >0 (5l) 
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and 



The other coefficients may he similarly evaluated. Once the coeffi- 
cients & 2 , . • • have been determined, all the perturbation- 

velocity components in the plane of the wing may be derived by means 
of the equations previously presented. 

A special case of equation (5l) gives equation (30) . The solu- 
tion for the swept trapezoidal wing may thus be obtained from the 
general equations that have been derived without the a priori assump- 
tion that the flow is conical. 

If only functions derivable from the x component of the pertur- 
bation velocity are required (such as pressure, lift, and drag coeffi- 
cients), equation (3), (6), (7), or (48) may be directly applied. The 
solution for the < wing of figure 4(f) that satisfies the Kutta- 
Joukowskl condition may be obtained by substituting equation (49) in 
equation (48) to yield 



For the special case of a wing with a straight leading edge (fig. 3), 
equations (53) and (8) give equation (12), which was previously 
obtained by intuition. 


RECAPITULATION 

A recapitulation of the formulas for the x component of the 
perturbation-velocity potential of a single thin flat-plate wing tip 


988 
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with arbitrary plan boundaries is given herein. The lift-distribution 
coefficient Cp^L of thin wing tips may be obtained from these 
formulas by the relation 

°P,L - (-2=p> - 1 § 

(a) In regions not influenced by the wing tip. 



(from equation 3(a) and fig. 1(a)). 

(b) In regions influenced by subsonic leading or trailing edges 
when no vortex sheet exists directly behind the subsonic trailing 
edge. 



(from equation (7a) or (48) and fig. 2(b)). 

(c) In regions influenced by subsonic trailing edges when an 
undetermined amount of vorticity is allowed in the vortex sheet 



(from equation (48a) and fig. 4(f)). 
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(d) In regions influenced by subsonic trailing edges when the 
Kutta-Joukowski condition is satisfied at the trailing edge. 


Sep Ua 
2(3* 


S- 


dv-du 


ab 


fv u )(v 7 ) 


(from, equation (53) and fig. 4(f)). 


Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, January 12, 1948. 
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APPENDIX A 
SYMBOLS 

The following- symbols are used in this report: 
a 0 , 8-\, ag, . . 


-p 

f(y) 

gO v ) 

k i 

k 2 

M 

m, n 

P 

S 

U 


u - § <«-*) 


v = 


_M 

2f 

M_ 

2P 


( I+Pt ) 


x, j, z 

I, n, t 

a 

0 

r 


a^.^ coefficients of power-series expansion of 
f(y) = &o + a l7 + a 2y 2 • • • 

pressure coefficient on top wing surface 

lift-distribution coefficient ( -2 C-q) 

perturbation-velocity potential on top surface 
of diaphragm in wake of so-called subsonic 
trailing edge 

undetermined function of v w 
constant greater than 0 

constant whose value lies between 0 and 1 
free-streara Mach number 
integer summation indices 

^ u 2" T w 

plan-form area 
free-stream velocity 


oblique coordinates whose axes lie parallel to 
Mach lines 


Cartesian coordinates 
angle of attack 

cotangent of free-stream Mach angle (ipci) 
gamma function 
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A 


a 


<P 


Subscripts : 

1 ) 2 y 3, • • • 

a 

b 

B 

D 

T 

w 

Examples : 

CTt 

*1 

CRD 


Sw(l+2) 

n 

r(n+i) 

cpa 


^a 


slope of streamlines near z = 0 plane measured 
in y = constant planes (diaphragm-surface 
slopes) 

slope of streamline (on wing slopes) with respect 
to z = 0 plane measured in y = constant 
planes 

perturbation-velocity potential on top wing 
surface or diaphragm 


numbered areas or curves 
portion associated with f (y) 
portion .not associated with f.(y) 
bottom (wing or diaphragm surface) 
diaphragm 

top (wing or diaphragm surface) 
wing 


slope on top wing surface 

slope of diaphragm in plan-area 1 

perturbation-velocity potential on top surface 
of diaphragm 

wing area 1 plus wing area 2 
curve v = vi(u) 
nl 

portion of perturbation-velocity potential on 
wing surface associated with function f(y) 

portion of diaphragm slope associated with 
function f(y) 


h 


value of | along leading-edge boundary curve 1 
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APPENDIX B 

DERIVATION OF EQUATION (6) 

The -velocity potential at point (x,y) for the wing tip of 
figure 2(a) is given by equation (5) (written in Cartesian coordi- 
nates) as 





The area a'b'd'c' may be subdivided into the areas fgkc', a'egf, 
and eb'd'k. The area fgkc' is the same as the area abed and 
the relation of the distances in the denominator Is the same for 
corresponding points. The first integral of equation (B2) may there- 
fore be written 
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(B3) 

On the other hand, the area b'od’ may be expressed as the 
area ghk plus the area edd'khg minus the area eb'd'k. The area 
ghk is the same as the area bod and the relation of the distances 
in the denominator is the same for corresponding points. The second 
integral of equation (B2) may therefore be written 



886 



If equation (Bla) is subtracted .from the sum of equations (B3) and (B4) and if expres- 
sions that lead to infinitesimals of order higher than the first in dx are neglected, the 
result is 



The first two integrals of equation' (B5) are surface integrals over the areas Sy^x 
and Sy g, respectively. Integration of the third and fourth integrals with respect to i over 
a strip of width dx removes one Integral sign and replaces d£ by dx. The width of the 
strip in the fifth integral along the £ direction may he determined with the aid of the 
following sketch of the area dd'k: 
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d| = 2(dx-dX2) 



where "the slope of the wing-boundary curve in the 

vicinity of the point d. Integration of— the fifth integral of 
equation (B5) with respect to | removes one integral sign and 
replaces &£ with the fourth member of equation (B6). Equa- 
tion (B5) thus becomes 


K 
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If equation (B7) is rewritten in oblique coordinates, equation (6) 
results: 
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APPENDIX C 

CALCULATION OF PRESSURE COEFFICIENT FOR WING 
OF FIGURE 3 BY METHOD OF REFERENCE 1 


When Ojj = -Op = ct, equation (5) gives for the velocity potential 
on the top surface of the wing shown in figure 3 


Now 


r^w 


cp = 2S 
^ Mrt 


du 


v w 


dv 


2XJa 

Mrt 

J 


u 2 (v w ) 

f Uw v v +k x u 

U2(v w ) 


'-k]U 




du 


2Ua 
. Mrt 


l^iV-u 

l\j [u w -U 2 (v v )] [v w +k 1 u 2 (v w )] 


(k^+v^) , k 1 [u w -u 2 (v w )] 

+ — .. — tan A . - ■ - ■ — 

v w +kiu 2 (v w ) 


^fkT 


1 


iMrflM?) 


Suhstitution of equation (Cl) into equation (C2) yields 


(Cl) 


(C2) 


Wo. J 

du2(v w ) 

k 

v w +kiu2(v w ) (k]_+l) , . 

ki[u v -u 2 (v w )J ^ 

0* 

1 " 3v w 

1 

uw-u 2 (v w ) ‘ ^ tan ^ 

v w +kiU2(v w ) 


(C3) 


« 
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The pressure coefficient on the top wing surface is then 

C = £aj 1 du2 ^ , Vw+kiu 2 (v v ) 

P Prt dv w } u w -u 2 (v v ) 

(04) 
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(d) 

Figure 4. - Continued. Integration boundaries for evaluating 
velocity potential in regions influenced by subsonic trailing 
edge 
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(f) 

Figure 4* - Concluded* Integration boundaries for evaluating 
velocity potential In regions influenced- by subsonic trailing 
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